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Abstract Theory of Random Matrix Ensembles have proven to be a useful tool in the 
study of the statistical distribution of energy or transmission levels of a wide 
variety of physical systems. We give an overview of certain ^-generalizations 
of the Random Matrix Ensembles, which were first introduced in connection 
with the statistical description of disordered quantum conductors. 

1. INTRODUCTION: RANDOM MATRIX 
ENSEMBLES 

With a few notable exceptions, the interaction between the community of 
mathematicians who work in special functions, in particular, those that are 
in the area of (/-series and basic Hypergeometric functions and the physics 
community has so far been minimal. In this review article, we will describe 
some developments in one area in physics, namely the Theory of Random 
Matrix Ensembles, where a (/-generalization has been proposed to be relevant 
for multifractal states near a critical regime. The term multifractal states 
refers to the novel fractal characteristics of the quantum mechanical wave 
functions of a disordered electronic conductor. It is well known in the physics 
literature that depending on the strength of the disorder, a quantum con- 
ductor may undergo a phase-transition from a metallic state which describes 
a good conductor to those of an insulating state where the (static) conduc- 
tivity in a finite sample is very small. In a region near the transition such 
multifractal wave functions are observed at least in numerical experiments. 
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In addition to giving an overview of the results obtained so far as well as 
providing some details that have not been published before, the hope of the 
authors is to draw the attention of mathematicians to some mathematically 
posed unsolved problems that are important from the physics point-of-view. 
We expect that this will also provide some additional insight into some of 
the asymptotic behavior of the (/-polynomials that may not be obvious from 
their conventional definitions. More precisely, we suspect that there exists 
certain "universal scaling limits," to be explained later, for at least a class of 
(/-polynomials in certain asymptotic regime that have not been fully explored 
by those mathematicians who are interested in (/-special functions. The work 
discussed here were done over last several years in many collaborations with 
Carsten Blecken, John Klauder, Vladimir Kravtsov and Nikos Nicopoulos. 

Originally, Random Matrix Ensembles (RMEs) were proposed by E. Wig- 
ner in the early sixties to describe the statistical properties of the eigenvalues 
and eigenfunctions of complex many-body quantum systems in which the 
Hamiltonians of such systems are only defined in a probabilistic setting [22]. 
Over the past fifteen years, they proved to be a useful tool in the studies of 
equilibrium and transport properties of disordered quantum systems, classi- 
cally chaotic systems with a few degrees of freedom, two-dimensional gravity, 
conformal field theory and chiral phase transitions in quantum chromody- 
namics [19]. The reason that such a wide variety of physical systems (with 
sizes ranging over more than nine orders of magnitude between a nucleus and 
a micron size disordered quantum conductor) can be described by theory of 
RMEs, which has no adjustable parameter in it, is that there exists a special 
double scaling limit in which the correlations between appropriately scaled 
eigenvalues in all these systems are universal, independent of any parameter 
that describes the particular system. 

In this section, we will briefly introduce the standard Gaussian RME, 
without any attempt to show how the expressions are derived. Detailed 
derivations can be found in [22] . The idea is to point out the reason for the 
existence of the zero-parameter universality under certain double scaling in 
Gaussian RME. We will then introduce the so-called "Coulomb gas" picture 
to just point out how the zero parameter universality had been "proven" 
to be valid beyond Gaussian RME. These will also serve the purpose of 
introducing notations and terminologies. In the following section we will 
show how this "proof" breaks down when the parameter q (0 < q < 1) is 
introduced. This then points towards to a new one-parameter universality. 
We will then review various known results as well as open problems for some 
of the (/-ensembles. Except for the section on circular ensembles, we will 
consider only unrestricted Hermitian matrices relevant for physical systems 
with broken time-reversal symmetry (which for example described a charged 
system in the presence of a magnetic field). The corresponding ensemble 
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is known as the unitary ensemble. We would like to point out here that 
although the Gaussian orthogonal and symplectic ensembles are by now 
very well understood largely due to the work of Tracy and Widom, the 
corresponding (/-generalizations have not been studied to date. 

1.1. GAUSSIAN RME 

Consider the set of all N x N Hermitian matrices M randomly chosen 
with the following probability measure 

P${M)dM = A N exp[- Tr(M 2 )]dM, (1.1) 

where An is a constant, Tr is the matrix trace and dM the Haar mea- 
sure. By going over to the eigenvalues-eigenvectors representation, it can be 
shown that the joint probability distribution of the eigenvalues X = (xi,i = 
1,2,... N) of the matrices can then be written (up to Vjy, the volume of the 
unitary group) in the form [22] 

N N 

P$(X)dx 1 ...dx N = Cg H {xi-XjfUe-^. (1.2) 

l<i<j<AT i=l 

Here the factor ni<j( x * ~ x j) 2 ar ises from the Jacobian of a change of vari- 
able, see [15]. 

The level correlations can be determined exactly by recognizing that the 
distribution can be written as a product of Vandermonde determinants of a 
set of (monic) orthogonal Hermite polynomials. The orthonormal Hermite 
polynomials satisfy the orthogonality relation 

oo 

J e~ x2 H n {x)H m {x)dx = 6 mn . (1.3) 



The main quantity of interest is the two-level kernel K^{x, y) from which all 
correlation functions can be obtained. For example the n-level correlation 
function is given by [22] 



N\ f f 

R n (x 1 ,x 2 , ■ ■ .x n ) = jjj— —y J ■■■ J P N (x 1 ,x 2 ,...x N )dx n+1 ...dx N 

— oo — oo 

= det[K N (xi, Xj)]ij = i, 2 ...n 



(1.4) 
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For the above Gaussian distribution the kernel has the explicit form 

N-l 

K%{x,y) = e -(* 2 +f 2 V 2 H n (x)H n (y) 

n=0 (1.5) 

-(x 2 +y 2 )/2 ntFJZ [H N (x)H N ^ 1 (y) - H N (y)H N ^ 1 (x)] 



VN/2 



(x - y) 



where we have used the Christoffel-Darboux formula for the last equality. 
In particular, the density of levels can be obtained simply as Kn(x, x). The 
Plancherel-Rotach asymptotics for the Hermite polynomials treat the large 
N behavior if x = y/2N + 1 cos 4>, or x = ^2N + 1 cosh tp, where 4> and oo 
are fixed and 

O<e<0<7r-e, < lo < ip < uj, 
see [34]. In the large N limit, with x and N as above, the density of levels 

a N ( x ) = Kg(x,x) ~ -y/2N -x 2 , (1.6) 

7T 

and zero otherwise. 

We will now consider the double scaling limit, given byN^oo,x,y^0 
such that the products x = (, y = rj remain finite. This par- 
ticular N dependence of the product is chosen so that the density a(Q = 
aN(x)dx/d( = 1 in the new scaled variable. We will ignore here rapidly 

falling density at the edges where it becomes zero and have focused 

our discussion to the middle of the density, i.e., x near where oat(O) = ^2K m 
In these new scaled variables, i.e., £ and rj and using the known asymptotics 
of Hermite polynomials, the two-level "limit kernel" becomes 

K G ((,v) = Sm[ (} C ~w *g(C)=^g(C,C) = 1- (1-7) 

7r(C - rj) 

This kernel, which is universal because it does not depend on N, gives rise to 
the famous Wigner-Dyson-Mehta correlations of the eigenvalues [22]. Note 
that if we choose the ensemble from a more general Gaussian distribution 

P N {M) = A N Exp[- Tr(aM 2 + bM + c)], (1.8) 

the parameters a, b, c can always be removed by a simple completing of 
squares. Therefore by an analogous rescaling the universality (parameter 
independence) of the two-level limit kernel is again reproduced. The details 
of this treatment are in Percy Deift's excellent monograph [15]. 
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1.2. BEYOND GAUSSIAN RME: THE 
COULOMB GAS PICTURE 

We now choose a non-Gaussian distribution of the ensemble defined by 

P%{M) = A V N Exp[-Tr(V(M))} (1.9) 

where V(x) is a suitably increasing function of x, so that the multiple inte- 
gral, 

oo oo 

J Pn(X) dx\... dxN, 

-oo 

exists for any fixed N, where the joint probability distribution of the eigen- 
values is 

N N 

P V(X) = C%H( Xi - Xj ) 2 l[e- V M. (1.10) 

i<j i=l 

By going through the same procedure as described above we can write the 
two-level kernel in terms of a set of polynomials that are orthogonal with 
respect to the weight function w{x) = exp(— V(x)) rather than e~ x which 
previously led to the Hermite polynomials. However, for an arbitrary V(x), 
one would not expect the corresponding orthogonal polynomials to have 
explicit representations, or simple properties, and therefore this is usually 
not a very useful description. Dyson introduced a Coulomb gas picture [16] 
by writing 

P%{X) = e? H , (1.11) 

and identifying 

N N 

H = - In \ Xi - Xj \ + - V(xi) (1.12) 

i<j i 

with the Hamiltonian of a "gas" of iV classical charges on a line interacting 
with each other with logarithmic Coulomb repulsion and held together by 
a confining potential V(x), at a "temperature" 1/(3. For our unitary en- 
sembles, (3 = 2. In the large iV limit, Dyson assumed that a continuum 
description is valid. (We will see later that this is where the (/-ensembles dif- 
fer in a subtle way) . This is a plausible approximation if the density of the 
charges increases with TV so that the spacing between the charges decreases 



/ 
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with N, and the N — > oo limit allows one to use a continuum approxima- 
tion. In this case, the two-level correlation function (related to the two-level 
kernel by eq. (1.4)) can be written as [5] 

^ / x 1 Sa(x) 

= ?Jvffi- (1 ' 13) 

The variational derivative can be evaluated in the coulomb gas picture by 
minimizing the energy (subject to the positivity of the density), which leads 
to the integral equation 

SV(x) + J Sa(y) In \x — y\ dy = constant, (1-14) 
J 

where J C R is the support of the density a(x). In general, the determination 
of J is a highly non-trivial problem and remain unsolved. The inverse of the 
equation then shows that the functional derivative of (1.13), and therefore all 
level correlations are independent of the choice of the form of the potential, 
so that the non-Gaussian ensemble has the same universality as the Gaussian 
ensemble provided that the assumption that J is a single interval is valid. 
This is the reason why many physical systems, which may be described 
by different constraints and hence different form for the distribution of the 
ensemble, still follow the universal zero-parameter description of the RME at 
the double scaling limit. Note that the proof depends only on the existence 
of the continuum limit where the electrostatic coulomb gas picture holds. 

1.3. BEYOND GAUSSIAN RME: THE 
ORTHOGONAL POLYNOMIALS 

We now reconsider the non-Gaussian case from the point-of-view of the 
orthogonal polynomials. As mentioned before, for any confining potential 
V(x), we can define the set of polynomials 4> n (x) orthogonal with respect to 
the weight function w(x) = e~ v ( x \ 



oo 

j , \ i . 



e v{x >(j) n (x)(p m (x)dx = 5 mn , (1.15) 



— oo 

such that the Christoffel-Darboux kernel becomes 

N-l 

Kl{x, y) = e -(vWv%))/2 J2 M*)My)- (1-16) 

n=0 

The point now is that although we don't know the polynomials for arbitrary 
V(x), the universality in the coulomb gas picture implies that polynomi- 
als orthogonal with respect to any weight function must have the scaling 
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property in the asymptotic region qualitatively similar to the Hermite poly- 
nomials, namely 

lim {-l) N N 1 l A H 2N {x) =7r- 1 / 2 cos7rC; 
lim (-l) N N 1 ^H 2N+1 {x) = Tr-^sinyrC, 

such that the double scaling that removes the ^-dependence from the density 
(to make it unity in the scaled variable) also removes the iV-dependence 
from the two-level kernel. In other words, the universality of RMEs for 
some confining potentials is related to some "universal" asymptotic scaling 
properties of the orthogonal polynomials. 

2. q-ORTHOGONAL POLYNOMIALS: FAILURE 
OF THE ZERO-PARAMETER 
UNIVERSALITY 

Clearly, various random matrix models can be defined either by choosing 
a confinement potential, or equivalently a set of orthogonal polynomials. 
It is therefore obvious to ask if the ^-polynomials share the same scaling 
characteristics as the classical polynomials. Historically, physicists working 
with random matrix ensembles did not know about (/-polynomials, so this 
question was never asked. It was a particular weight function arising in 
the problem of disordered conductors [30,32,33] that eventually led the first 
and second named authors to seek help from Prof. George Andrews who 
suggested that Mourad Ismail be brought on board, starting a very fruitful 
collaboration. 

In the next section, we will consider some specific examples of the q- 
ensembles. Here we mention the important general features of g-polynomials 
that are relevant for the breakdown of the double scaling universality for the 
g-ensembles. 

For classical orthogonal polynomials, the density <jn{x = 0) scales with 
some power of N (e.g., a^{x = 0) oc y/N for Hermite polynomials). It has 
two consequences. First, the spacing between nearest two levels scales as a 
power of 1/-/V, so that the continuum limit can be achieved by taking the 
N — ► oo limit. Second, the fact that lim o>/(0) — > oo guarantees that if 

the moments of the weight function 



oo 




— oo 



are given then w is uniquely determined. For those g-polynomials that cor- 
respond to the indeterminate moments problem, the density crjv(O) is in- 
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dependent of N. The fact that in this case lim <tjv(0) is finite implies 

N^oo 

that the corresponding moment problem is indeterminate [1]. Clearly we 
expect the asymptotics of these polynomials to have different scaling prop- 
erties. We have observed however that all the cases of g-orthogonal poly- 
nomials when the measure (weight) of orthogonality is not unique one can 
always find a weight function w(x) so that (\nw(x))/ In 2 x has a limit as 
\x\ — ► oo. This probably contains a new universality law. This distinguishes 
g-orthogonal polynomials, that is polynomials whose recursion coefficients in 
the orthonormal form grow exponentially, from other cases of indeterminacy 
with recursion coefficients growing at a polynomial rate. 

It turns out that the weight functions for some g-polynomials behave as 
g-cin (x) £ or i ar g e x _> qo f or some positive constant c. An indication of this 
can be seen from the growth of the the recurrence coefficients as in the con- 
jecture formulated at the end of this section. If the recurrence coefficients of 
a symmetric set of orthogonal polynomials grow as 0(q~ aN ), for a > 
then the V(x) = 0((lnx) 2 ), see [20] for an example. As the moment problem 
is indeterminate, there are different choices of the weight function possible 
for the same recurrence relation. However, V(x) oc ln 2 (x) for large x is a 
characteristic property for all choices. Comparing with the Coulomb gas 
picture, this means that the confinement potential in this case is "weak," 
comparable to the repulsive interaction ln|x — y\, as opposed to a "strong" 
power law confinement in case of classical polynomials which always domi- 
nate over the logarithmic repulsion. Note here that in order for the moment 
problem to be determinate V(x) should grow at least as fast as | x | as x 
tends to infinity. This difference is crucial. In fact, if we choose a weak con- 
finement potential of the form V(x) = ln 2 (l + ax) for which no explicit (or 
closed) form of the orthogonal polynomial is known, the level correlations 
obtained by numerical methods turn out to be the same as those obtained 
analytically from exact solutions of the g-models [11]. Thus the soft con- 
finement, N- independent spacing of the levels and indeterminate moments 
are all related to each other, and are generic characteristics of what we will 
call the "(/-Random Matrix Ensembles." We will see later that these ensem- 
bles are related to the physics of multifractality near the critical regime of a 
phase transition where the wave function of a particle in a quantum conduc- 
tor changes from an extended state to an exponentially localized state when 
randomness is increased. 

The following conjecture is based on observations we made over the years. 



Conjecture 1 Let p(x) be symmetric orthonormal polynomials satisfying 

(X) 
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and assume that a n q an has a limit as n — > oo, for some a > 0. If w is a 
weight function for the p n 's then 

\nw(x) ~ — cln 2 |x|, for some c > 0, 

as x — > oo or x — > — oo, frtii not necessarily for both. 

In fact we even believe that the conjecture holds even for discrete measures 
in the sense that the w(x n ), the mass at x n satisfies the asymptotic relation 
lnw(x n ) ps —cln 2 \x n \. There are ample examples to illustrate this behavior, 
sec [20]. 

2.1. THE q-HERMITE RME 

Following [27], we will define the q-Hermite random matrix model by the 
q-Hermite polynomials orthogonal with respect to the Askey weight function 
w A [3] 



h n (sinhu; q)h m (smh.u; q)wA(u; q) du 

= q -n(n + l)/2^ qUq . q)oQ Hl/q)6n 



(2.1) 



where 



WAM = (- q e-*- q e*»; q )J ° <q<1 ' M 



Here — oo < u < oo, and 

n n 

(a;q)n = II (1 - aq 1 " 1 ) ; (oi, a 2 , . . . a n ; q) m = JJ(aj, q) m , (2.3) 
i=i i=i 

where m = 0, 1, . . . or oo. Ismail and Masson [20] referred to the polynomials 
h n as the q _1 -Hermite polynomials. Note first of all that the polynomials in 
(2.1) are in the variable sinhn, which means that the confinement potential 
has to be compared in the variable x = sinh-u. To check the asymptotic 
behavior of the potential, we use the Jacobi triple product formula [18] 

oo 

(q, -zy/q, -y/q/z, q)^ = £ q n '' 2 z n . (2.4) 



Choose z = q 1 l 2 e 2u , giving 
1 + e~ 2u 



w A (u;q) (q,q) c 



^ g n(n+l)/2 e 2n_ ^ 
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Note that the right hand side of (2.5) is a constant multiple of a theta 
function. Also bear in mind that as a function of x the weight function 
should be wa (u;q)/ cosh u, and x = sinh u. So applying the Jacobi imaginary 
transformation, we find after going back to the x variable, 

V A (x-q) =0((sinh- 1 x) 2 /7), (2-6) 

which is O(lnx) 2 . Here 7 = ln(l/g). 

From the asymptotic formulas derived by Ismail and Masson [20], the 
Christoffel-Darboux kernel is found in the limit N — ► 00 to be 



y/w{u)w{v) 

K A {u,v;q) = 



[(g;g)oo] 2 ln(l/?) (2.7) 
x {-qe u+v ; -qe'^; qe u ~ v ; -qe~ u+v ; q) ^ . 

Note that the density, in this ./V — > 00 limit, is 

<j A (u;g) =i^(u,u ;(? ) = l/ln(l/g). (2.8) 

We note here that the Christoffel-Darboux kernel has a limit as N — ► 00 for 
fixed u and i> and therefore the density also has a limit. Compare the above 
with the ordinary Hermite case 



a N (u = 0,q= 1) = . (2.9) 

7T 

Note that for any q < 1, the iV — > 00 limit of the kernel exists. This is a 
discontinuous change, and does not allow us to study the important question 
of how the zero-parameter universality breaks down in a physical system. 

It is tempting to conjecture that the physically appropriate generalization 
of the density from the q = 1 case should involve 

W « S T^ = ^^ ; ' se " T (2 - 10 > 

replacing N, which is equal to iV for 7^ <C 1, but becomes independent 
of N for 7 AT 3> 1. This means that we should explore the asymptotic 
behavior of ^-polynomials in the double scaling limit where 7 ^ 0, iV ^ 
00, but the product remains finite. Although they have been used in 
certain special cases [8], systematic studies of such double scaling limits 
of (/-polynomials have not been done. The kind of limits just described; 
q becoming iV dependent would require analysis involving varying weight. 
This should be of interest to mathematicians who work in the asymptotics 
of orthogonal polynomials. 

The kernel (2.7) is independent of N, but depends on q. The question 
then becomes, can the g-dependence be scaled out by changing to variables 
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such that the density becomes unity, thereby restoring the zero-parameter 
universality? The answer was given in [27]. Here we provide some detailed 
steps of the calculation (and supply the g-dependent constant factors). We 
first express the infinite products in terms of Jacobi theta functions B\ [18], 

V ' Joe 2(q;q) OQ ismh((u — v)/2) 

x (-qe u+v ;-qe- iu+v) ;q) (2.11) 

1 e 2 (i(u + v)/2;^g) 
2q 1 / s {q;q)oo i cosh((w + v)/2) ' 

The kernel can be expressed as 

K( Q{u,v)Q 2 (u,v;^q) 9i{i{u - v)/2; y/q) 

KA[u,v,q) - %!)oo]Sh(1/!)?V8 iB mh((«-t,)/2) ^ 

where 



Vcoshucoshu 9 2 (i(« + u)/2;v/9) , . 

S2(«, v) = — ri—'i ^>2(u, v; Jq) = — . (2.13) 

1 ; cosh(^ + V )/2' 21 ' Vyj V^KVP^^ 1 ^ 

We now use Jacobi's imaginary transformations for the theta functions [36], 

V 1_ V 7 ' (2.14) 

where p = e ~ 271-2 / 7 . The kernel now becomes 

<9 4 (7r(ii + v)/t,p) 9i(ir(u-v)/T,P) 



K A (u,v;q) = a(q)Q(u,v) 



where 



v / 6»4(27rn/7;p)6» 4 (27r?;/7;p) sinh((« - v)/2) ' 

(2.15) 



a(q) = \[—-—— u q \\^ . (2.16) 

W V 7 2 7 [fe<?)oc] 3 1 ^ 

Defining new variables 

C = Vt; f] = v h ( 2 - 17 ) 
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we finally have the appropriately scaled two-level kernel 



Ka{(, V, q) = a(q)n(jC, 7r?)04(C, v) 7^A 



XC - v);p) 



sinh((C - 77)7/2) ' (2.18) 

if A (C,C; <?) = i, 



where 



e 4 ((,) = T M^lL (2.19) 

Apart from the prefactor a(q), this is the limit kernel quoted in [27]. Note 
that the kernel remains a function of 7 = ln(l/g), so that the distribution 
is no longer parameter independent. However, the dependence on the levels 
as well as q are quite complicated, and various limiting cases turn out to be 
quite interesting. In the rest of the section we will consider some of these 
limits. 

One corollary of the above expressions for the kernel before and after the 
imaginary Jacobi transformations is that we can equate the density Ka(u, u) 
before and after the transformations to obtain the relation 

[(q-qUf = (^yVV/ 4 [(PV)J 3 . (2.20) 

Using this, and the product expansion of the Jacobi theta function 

00 

e l {u;p)=2(p 2 ; P 2 ) ooP 1 / 4 smuYl[l-2p 2n cos2u+p 4n ] , (2.21) 

1 

we can now consider some characteristics of the limit kernel. 
For £ = 77 + e, we get 



— 7 sin[7rej i r 

KA{e) ~ 27 sinh[ 7 e/2] U 



V" .,.0 



71=1 



1 + -, ^ — ttt sin me] 

(l-p 2n ) 2 1 1 



+ o(l) (2.22) 



This expression depends only on the difference ( — 77. We will see that the 
lack of this "translational invariance" for the general expression (2.18) has 
important consequences. For p <C 1, i.e., 7/271" 2 <C 1, the second term in the 
product can be neglected, and we get the simple hyperbolic kernel 

^-„)^ si :;^ 2) + o (1); (, 23) 

The above limit kernel for the g-Hermite model for j/2tt 2 <C 1 has been 
shown to be equivalent [21] to a banded random matrix model that describes 
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multifractal wave functions at the critical regime for a quasi one dimensional 
disordered quantum conductor with weak multifractality [23]. Comparing 
the two, the parameter 7 = ln(l/(/) can be explicitly identified with the fractal 
dimensionality at the critical regime [21]. In contrast, the universal kernel 
resulting from the ordinary Hermite polynomial describes chaotic systems 
with extended wave functions, i.e., those wave functions that extends to 
the edge of the sample. Thus the (/-deformation is in some way related to 
the physics of multifractality, although the nature of the relation is not at 
all clear. One major problem in understanding the connection is that we 
do not know how to study the crossover from q < 1 to q — > 1~. For a 
fixed q < 1, the large N limit kernel is independent of N. This prevents 
us from exploring in details how the double scaling fails as one approaches 
the limit q — > 1~ when the (/-Hermite polynomials reduce to the ordinary 
Hermite polynomials and the restoration of the zero-parameter universality. 
As pointed out before, an interesting direction is to study the double scaling 
regime 7 — ► 0, N — ► 00, keeping ^yN finite, and then consider the limit 
kernel. Properties of (/-polynomials in such special double scaling limits 
have not been explored yet. 

The situation is completely different in the case where 77 is near — £. The 
physical importance of this "off-diagonal" correlation was emphasized in [11]. 
Consider i] = Again for p <C 1, the ^-functions simplify, and we obtain 



Since coth7^ — > 1 for large £, the kernel have long range oscillations of 
frequency 2n whose amplitude grows linearly with 7 (as long as 7 <C 2tt 2 ). 
The case for larger 7 should still have the oscillations, although the prefactor 
will change. In other words, the kernel exhibits long range correlations in 
the 77 = — £ direction. Another physically interesting limiting case is the 
kernel in the limit q — > 0, i.e., 7 — ► 00. The kernel was heuristically shown 
to be of the form [9, 21] 



in squares of unit size along the two diagonals ( = n and ( = —n, where [2(] 
is the integer part of 2Q. For all other values of £ and 77, K(£, rf) = 0. This is 
a surprising result. It is not clear if this result is valid only for the q- Hermite 
polynomials, or is it a generic feature of (/-polynomials in the q — > limit. In 
particular, the asymptotic properties of the (/-polynomials for large N in the 
scaling limit N — > 00, 7 — > 00, keeping the ratio finite would be interesting 
from physical point-of-view. 



Ka(C, -0 = 7T coth( 7 C) sin(27rC) + o(l). 



(2.24) 




(2.25) 
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2.2. AN ALTERNATIVE g-HERMITE RANDOM 
MATRIX MODEL 

Since the moment problem associated with the g-Hermite polynomials is 
indeterminate, it is instructive to consider a different g-Hermite model to 
check if the above results are sensitive to the particular choice of the weight 
function. One such model is to choose 

oo 

J h n (x; q)h m (x; q)w qH {x; q)dx = (2.26) 

— oo 

where w q n is the weight function onsidered in [4] 

w qH (x;q) = e- 2 ^ h ~ lx ^^. (2.27) 

This model of g-RME was first considered in [26]. Here we provide some 
details. 

Note that the asymptotic behavior of the confinement potential in this 
model is the same as for the Askey weight. Equation (2.26) can be rewritten 
in the variable ( = sinh -1 ar/7 



/i n (sinh7C; q)h m (smhj(; q) cosh^e' 2 ^ 2 d( 

(2.28) 

~ (q,q) n x 



2 7 g (n+l/2)V2 " nm ' 

to take advantage of the Ismail-Masson result for the kernel, which becomes 

= 7 ^ V^oshjC cosh 7?7e -7(C 2 +» 2 ) (2.29) 

V vr (g,g)oo 

X (-qe^+rt ; - qe -i«+v) ■ qe^~^ ; qe'^'^ q) c 



loo- 



In terms of Jacobi theta functions, the kernel has the form, after taking the 
N — > 00 limit, 

— fin 1 

K qH {(,v;q) = \ — 



vr 2iq 1 / 8 {q,q) 00 

V 7 cosh 7C cosh 777 (2.30) 



cosh7((" + r/)/2sinh7(C — rj)/2 
x e 2 MC + ri)/2; v^ei[i 7 (C " v)/2\ Vq}- 
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Using the Jacobi imaginary transformation as before, we obtain 

KMC m q) = c(g)n(7C, vi)o*{*{< + v)) 2^-^/2] ' (2 ' 31) 

where 

= vf W^P (2 ' 32) 

and fi has been defined before. Note that K q H((, () = #4(27rC) is not unity, 
so one needs to rescale the variables again. However, for 7 <C 2-7T 2 , the density 
is approximately unity, so we can compare this kernel with the one from the 
Askey weight, which is the same in this regime. Moreover, it is also the same 
for C = —f]. Heuristic arguments shows that the q — > behavior to be similar 
for all potentials that have the same large x behavior [9], so we expect it 
to be the same as well. Note that in terms of the original variable x, the 
confining potential had an oscillatory component in the Askey weight while 
the density was monotonic; in g-Hermite weight the potential is monotonic 
while the density is oscillatory. Thus, the two models are very different. 
It is only after the appropriate scaling that the kernel, and therefore the 
physical properties, show similar behavior. Note on the other hand that the 
(/-Hermite weight is related to the Askey weight as 



/ v. 27^/2(1 9)00 
w A (x; q)dx = — e A (2nQ-p) w i H ^ Q)d(- (2.33) 

Since the Jacobi theta function is periodic with respect to a shift #4(2:; q) = 
0^(z + 7r; q), the two weights differ only by a periodic factor with respect to 
the shift £ — > £ + m/2 for integer m. Whether the physical correlations are 
always equivalent for all choices of the weight function that correspond to 
the g-Hermite polynomials, is an interesting open question. 

2.3. THE q-LAGUERRE RME 

We have already seen that the appropriately scaled physical correlations 
are independent of the choice of the weight function within the g-Hermite 
random matrix model and appropriate parameter regimes. In order to see 
if these are generic properties valid for other ^-random matrices, we now 
consider as an example the g-Laguerre RME, first considered in [12,13]. It 
is defined by the weight function [24], 
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with the orthogonality relation 



oc 

/ 



L n (x; q)L m (x; q)w qL (x; q)dx = ^^^ n Km, (2.35) 



where L n (0;q) = 1. The orthogonality is on x > 0, so the boundary at 
x = makes the model qualitatively different from the g-Hermite model. 
Moreover, as shown in [12,13], for small x the confining potential is linear in 
x as opposed to quadratic dependence for the g-Hermite potential. However, 
the large x behavior is again In 2 x, so the breakdown of the zero-parameter 
universality should again be similar. In particular, note that the density at 
x = can be obtained directly for any N, 

1 _ p-iN 

a% L (x = 0;q) = (1 - q N )/ln(l/q) = . (2.36) 

7 

Compare this with (tn{x = 0; q = 1) = N. As we suggested for the g-Hermite 
model before, studying the double scaling limit of the g-Laguerre polynomials 
7 — > 0, N — ► oo with the product "yN finite, would allow us to study the 
crossover regime from the iV-dependent density and the TV-independent q- 
model. The g-Laguerre two-level kernel was first written down in [7], and 
here we will provide some details. We use the asymptotic properties of the 
g-Laguerre polynomials [12, 13] 

N+l 

L N (x; q) = L^x; q) + (a; q) - L^x/q; q)} + 0(q N+2 ) (2.37) 
where 

Loo(x; q) = 4 2) (2y/x(l - q); q) . (2.38) 

Here jjp(z;q) is the g-Bessel function. Note that as before, the N — ► oo 
results are independent of N . The kernel requires the combination 

L = L N+ i(y; q)L N (x; q) - L N (y; q)L N+1 (x; q) 

= q N+1 [Loo {y/q; q) L^x; q) - L^y; q)^ (x/q; q)} (2.39) 
+ 0(^ +1 ). 

Using the asymptotic expression for the g-Bessel function for q <C 1, [14] 

4 2) (2y/x(l-q); q) ~ ^cos {(tt/t) In y^} 

2 77- 2 ( 2 - 4 °) 



x exp{[hr y/Qx - tt 2 /4\/^} 
Q:=l-q, 
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we get after some algebra 



L = 4?F ^Q c -^ 2 /7+7/4 c [ln 2 Qx+hi 2 Qy] 



1 



x \_-\fx cos {(71-/27) In Qy} sin {(71-/27) m Q x } 

- sin { (7r/2 7 ) In Qy} cos { (tt/2 7 ) In Qx}} + o(l) . 



(2.41) 



For large x such that q <C Qx, we can use Jacobi's triple product identity 
above to obtain 



—00 v ^ 



(2.42) 



Approximating the sum by an integral we get 




[w qL {x; q)} = \j -j-e ^ [In Qx + 7 /2] 2 + o(l). 



(2.43) 



Using these expressions, we obtain the limit kernel 



K qL {x,y) 



6(9) 



x-y 

•y\l/4 



x\ 1/4 
— sin 

yj 



7T 
L27 



In 



cos 



7T 
27 



InQy 



sin 



7T 
27 



In 



cos 



7T 
27 



InQx 



(2.44) 



where 



6(g) = 



8ir ( p 



1 



3/2 



2\ 1/8 



(2.45) 



and < 7 < 2ir. Introducing the variables ( = ^^-InQx and r) = 
In Qy, the kernel becomes 

7 sin(C/b) cos(7//b) - e"^^ sin(??/b) cos(C/b) , , 

7r smh[(C — 7?)7/7r6J 



and K q L((,() = 1 + ( 7 /27r) sin(2£/6). For 7 <C 2ir, the density is approxi- 
mately unity and independent of ( and we can compare the kernel with the 
q-Hermite model. For £ ps 77, K q L(Q,rj) has the same form as ^^(C) 7 ?) f° r 
the g-Hermite model, and therefore the level correlations are similar. 
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2.4. 



q-CIRCULAR RME 



For mathematical simplicity, Dyson introduced the Circular Ensemble 
where the eigenvalues lie uniformly on the complex unit circle [17]. A pos- 
sible ^-generalization appropriate for eigenvalues of scattering matrices for 
disordered systems was introduced in [28] defined by the weight function 

2 



w s (9;q) 



1 

27 



(q 1/2 e ie ;q) c 



(aqV*e*;q) c 



< q < 1, a q < 1. 



(2.47) 



This defines the Szego polynomials generalized by Askey [34], orthogonal on 
the unit circle 



2tt 



1 

2^ 



^ m (e i6 )^ n (e i6 )w(9;q)d9 = 5 r , 



(2.48) 



where the overline denotes complex conjugate. The two-level kernel is then 
given by 



N-l 



K S (0, <P)Vws(0; q) vW^ Q) E ^(e ie ; q) <f>n((e^ ;<? ). (2.49) 



n=0 



In the large N limit, the density was obtained in 



N 

<0;q)^- 



1 + 



1 



(1 - q)N 



In 



1 — la^fq cos 6 + a 2 q 
1 — la^fq cos 9 + q 



(2.50) 



There are two important features of this result. First, the density has finite 
N correction to the uniform density N/2ir of the a = 1 Dyson circular 
ensemble only in the q — > 1 limit such that the product (1 — q)N is kept finite. 
Thus the physically interesting limit is again a double scaling limit, although 
different from the previous ones. The second point is that although properly 
normalized, the density becomes negative for sufficiently small values of 9 < 
9 C , where 9 C depends on the parameter a. (Physically, more disorder means 
larger 9 C ). It was pointed out in [28] that this negative density arises as 
a result of the opening up of a gap, which requires that this possibility be 
included from the beginning by demanding orthogonality of the polynomials 
on the appropriate range 



C 



2tt- 



^ / * m (e i0 )* n (e*)w(e;q)d6 = S n 



(2.51) 



These are no longer the generalized Szego polynomials. The problem sug- 
gests considering orthogonal (/-polynomials defined on parts of a circle. 
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2.5. BI-ORTHOGONAL g-RME 

The bi-orthogonal random matrix models were proposed in [25], where 
the joint probability distribution is of the form 

N N 

P N (X) = C N Yl [r (x m ) - r (x n )} [s (x m ) - s (x n )} ]J w (x n ) . (2.52) 

m<n n=l 

Here r(x) is a polynomial of degree h and s(x) is a polynomial of degree k. 
The kernel is then given by 

N-l 

K$(x,y) = y^(x) y/My) £ Y n (x)Z n (y) (2.53) 

n=0 

where Y and Z are polynomials in r(x) and s(x) respectively, satisfying the 
biorthogonality relation 

J Y n (x)Z m (x)w{x)dx = g n S mn . (2.54) 
I 

Here g n is the normalization constant, and / refers to the appropriate range. 
It reduces to the ordinary ensembles for r{x) = s{x) = x. For the special case 
of r(x) = x, s(x) = x k , the biorthogonal limit kernel for the ordinary Hermite 
and Laguerre ensembles were recently obtained by Borodin [10] in terms 
of Wright's generalized Bessel function. For the special case of r(x) = x, 
s(x) = x k and w(x) = w g L(x;q) (the g-Laguerre weight defined earlier), the 
polynomials are known as the g-Konhauser biorthogonal polynomials [2]. 
However, the asymptotic form of the limit kernel for this model for large N 
originally proposed in [25] has not been obtained yet. 

3. OTHER PHYSICAL PROPERTIES 

While the two-level limit kernel in principle determines all physical cor- 
relation functions for unitary ensembles, evaluating a physically interesting 
property is not always straightforward. We mention two such properties 
considered for g-RMEs so far. 

3.1. FREDHOLM DETERMINANTS FOR 
q-RMES 

An important quantity for studying the spacing distribution of adjacent 
levels is the Fredholm determinant associated with the limit kernel. Tracy 
and Widom [35] considered Fredholm determinants for general random ma- 
trix models described by orthogonal polynomials. They showed that under 
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quite general conditions, one can obtain a system of partial differential equa- 
tions for ensembles defined by a weight function w(x) if x k w(x) is bounded 
for k = 0, 1, . . . and w(x) is continuously differentiable within the interval. 
Nishigaki [31] has used this approach to obtain the diagonal resolvent of 
the limit kernel (2.23). This is valid for 7 <C 2tt 2 , and has been shown to 
describe the spacing distribution of disordered conductors near the weakly 
multifractal metal to insulator transition region. (An earlier attempt to de- 
scribe the spacing distribution was also limited to the weak multifractality 
regime [29].) The corresponding result for the more general case valid for all 
q has not been obtained yet. Numerical results show a novel type of limiting 
"triangular" distribution in the limit q — > [9]. 

3.2. PARAMETRIC CORRELATIONS 

In the presence of an external field X, the "parametric" level density 
correlation function for g-Hermite model was introduced in [8] 

S (x, y; X, X') = cw (y/cx) w (y/cy) 

00 N-l 

x ^2 ^2 K (sinh (y/cx) ; q) h n (sinh (y/cy) ; q) 

n=N m=0 

x h rn (sinh {yTcx) ■ q) h m (sinh (y/ty) ; q) ^-n^X-X' f 

(3.1) 

where c = ^ir 2 /2N q and to = c/fi is related to a disorder parameter. In the 
limit X = 0, this reduces to the ordinary two- level correlation function 

S(x, y,X = 0) = R 2 (x, y) = a{x)a(y) - K\x, y). (3.2) 

In the limit q = 1 this reduces to the usual parametric density correlation 
function [5]. The expression (3.1) was considered in the special case x = y = 
only, in the double scaling limit 7 — ► 0, iV — > 00, with the product 7^ 
finite. The general case for all x, y is not known. Other physically important 
quantities are the "velocity correlator" defined as 

C(E l ,E 2 ;X,X') =JI ^L^S(x,y;X,X') dxdy, (3.3) 

and the 'parametric number variance' defined as 

U (E 1 ,E 2 ; X) = J J [S(x, y; 0) - S(x, y; X)} dxdy (3.4) 

Again only the limit E\ = E 2 = has been considered in [6,8]. 
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4. SUMMARY 

The two-level kernels (and other correlation functions) for (/-ensembles 
require the asymptotic behavior of various orthogonal (/-polynomials and 
their weight functions for finite N in various double scaling limits. Some of 
these scaling properties seem to share common features in certain range of 
parameters for apparently very distinct (/-polynomials. 
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